First Passage Time Statistics For Systems Driven by Long Range Gaussian Noises 
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We examine the mean first passage time for a particle driven by liiglily correlated Gaussian 
fluctuations to reach one or more predetermined boundaries. We discuss a numerical algorithm to 
generate power-law correlated fluctuations and apply these to three physical examples. One is the 
• arrival of a free particle at either end of an interval. The second is the decay of a particle from 

an unstable state. The third is the time for a particle to cross a barrier separating one well from 
another in a double well potential. In each case a comparison with the first passage time for a 
particle driven by Gaussian white noise is presented, as is an analysis of the dependence of the first 
passage time properties on the correlated noise parameters. 



I. INTRODUCTION 



The arrival of a random process at a particular state often triggers some important behavior. Among the large 
number of instances that one can name are the firing of a neuron, the nucleation of a phase associated with a phase 
transition, the triggering of an alarm, the occurrence of a major earthquake, and the crossing of an activation barrier 
by a reaction coordinate that converts reactants to products in a chemical reaction. These examples range from the 
^ ' macroscopic to the microscopic, indicating that the arrival problem is of interest on all scales. In many instances it 
is not only important to know the statistics of occurrence of such events, but more specifically the statistics of first 
occurrence. This then leads to the study of the statistical properties of the time that it takes a random process to 
^ ] reach a specified state for the first time, that is, the mean first passage time. The first passage time literature is 

• enormous and extends over many decades ||, |[ |[ |j . 
00 ] For simplicity we restrict our attention in this work to one-dimensional random processes x{t), although the gen- 
. eralization to vector processes or even to fields is conceptually straightforward. The statistical properties of the 

■ temporal evolution of the random process are typically described in one of two ways. One is by way of an evolution 
equation for the probability distribution P{x, t) for the random variable to take on the value x{t) = a: at time t (the 
random variable and the values that it can take on are often denoted by the same symbol, x in this case, although 
the distinction should be kept in mind). The other is by way of an evolution equation ("Langevin equation") for 

' the random variable itself. From the solution of the appropriate evolution equation (via steps that are in general 
non-trivial) one can then obtain the first passage time statistics for x{t) to reach a prescribed value for the first time. 
I In either case, one often thinks of the evolution of x{t) as being driven by fluctuations 7]{t) of prescribed statistical 

■ properties from which the statistics of x{t) then follow. 
Ch , Let us then turn to the "driving fluctuations" r]{t) and their typical statistical properties. By far the most common 

assumption is that fluctuations have a Gaussian distribution, although there has always been a great deal of interest in 
fluctuations that are distributed in some other fashion. The most common examples include Levy distributions |^ 
(distributions with very long tails), and dichotomous processes ^ |^ (where the random variable can take on only 
two values). The former typically arise when the fluctuations are a result of many multiplicative inputs; the latter 
often serve as an analytically tractable prototype. 

In order to specify the fluctuations fully, one must also explicitly state their correlation properties. If they are 
Gaussian, 

P(7?) = ^e-'V2.^^ (1) 



then only the average {ri{t)) (usually and here as well taken to be zero) and two-time correlation function (assumed 
stationary) j{t — t') = (r]{t)ri{t')) need to be specified, since all other correlation properties then follow. In Eq. (|l|) 
(T^ = 7(0). For other distributions it is in general necessary to also specify higher order correlations functions, even 
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though these are often difficult to determine. The early literature dealt primarily with 5- correlated or "white noise" 
processes, 

7(i) = 2DS{t) (2) 

and D (the integral of the correlation function) measures the intensity of the noise. In this case the second moment 
tJ^ diverges as D/At, where Ai ^ is the "width" of the (5-function and therefore the shortest time scale in the 
problem. The evolution equation P(x, t) for a random process driven by Gaussian (5-correlated fluctuations rjU) is the 
familiar Fokker-Planck equation ||l|, , and the associated first passage time properties are well understood [Q, ^ ^ . 

About two decades ago a great deal of attention began to be directed toward understanding the effects of "colored 
noise" , that is, of driving fluctuations that are not (5-correlated. Most of the attention was focused on exponentially 
correlated noise, 

lit) = (3) 

T 

where r is the correlation time of the noise. The associated evolution equation for P{x, t) and related first passage 
time properties as a function of the correlation time were studied in detail and are also well understood For 
instance, it is firmly established that the mean first passage time from one potential well to another of a process x{t) 
driven by exponentially correlated noise increases with increasing correlation time. It is also understood that the 
particular exponential form of the correlation function is not crucial in the qualitative features of the first passage 
time statistics. The most important feature determining the qualitative features of the first passage time statistics is 
the fact that there is a finite correlation time associated with the fluctuations. 

More recently there has been considerable interest in (typically Gaussian) fluctuations that display long-range 
power-law correlations with an infinite correlation time. Such highly correlated fluctuations have been considered in a 
broad array of circumstances ranging from the biolo gica l to the physical, from the economical to the atmospheric, and 
encompassing theoretical and experimental studies |l^ . There are two separate issues that need to be addressed 
with such highly correlated noise. One is the numerical generation of the correlated noise itself. Only recently have 
the traditional Fourier filtering methods been superseded by far more efficient procedures |l^, |lj, |l^, |l6|, that will 
briefly be reviewed below in the context of our systems. The second issue is that of the effect of such a driving noise 
on the first passage time properties of the system that is driven by these highly correlated fluctuations, and this is 
our principal interest in this paper. Herein we examine the mean first passage time of a particle, driven by highly 
correlated fluctuations, in three different physical situations. In the first case, we discuss the arrival of a free particle 
(moving superdiffusively) at an absorbing boundary. As a second case, we study the decay of a particle from an initial 
unstable state. Finally, we discuss the time for a particle to cross a barrier separating one well from the other in a 
double well potential 
In SecJ] 
0^ 



IV, and 



we summarize the main steps in the generation of long ranged correlated Gaussian noises. In Sees. HI 



we present respectively the arrival of the superdiffusive particle at an absorbing boundary, the decay of an 



unstable state, and the barrier crossing problem. Finally, we conclude with a summary in Sec. VI 



II. GENERATION OF LONG RANGE CORRELATED GAUSSIAN NOISE 

The main steps to generate noises with a Gaussian distribution and arbitrary correlation properties can be found 



in Refs. |12, 14, hSL 19, [L7[. Here we present a summary of these "spectral methods" including some details relevant 



to our applications. The goal is to generate a Gaussian noise ri{t), with correlation function 7(t) defined by: 

{^{tUt'))=^{t-t'), (4) 

and with a Fourier transform 



7(0.)= / dte-^-'jit). (5) 

J —oo 

This correlation function may be specified analytically or numerically. In the lu Fourier-space, the transformed noise 
rj{uj) has a correlation function 

{?j(uj)?j{uj')) = 2TT-f{Lj)S{uj + oj'). (6) 

The algorithm for the noise generation can be summarized as follows. First, the time interval {0,t) is discretized 
into N = 2"' intervals of mesh size At. This time interval has to be much smaller than any other characteristic time 
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FIG. 1: Comparison of Eq. (|l^) (circles) with correlation functions generated with the algorithm described in the text (solid 
lines). Upper curve: /3 = 0.3 with to = 0.0026. Lower curve: /3 = 0.6 with to = 0.0022. Other parameters for both curves: 
e = 20.0, N = 2^'-*, At = 0.01, loq = 0.0001. These simulation resuhs are averaged over 100 realizations. 
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FIG. 2: Same as Fig. hibut on a logarithmic scale that highlights the power law behavior. 



of the system (and hence N must be sufficiently large) because it is used as the time integration step in the numerical 
simulation of the stochastic evolution equations. The intervals in time will be denoted by a roman index and the 
resulting frequency intervals in Fourier space are denoted by a greek index. In discrete Fourier space the noise has 
a correlation function given by 



{'q{uj^)ri{uj^>)) = 7VAi7(w^)(5^+^/,o, 



(7) 
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where ^{ujf^) can be constructed as follows: 



t]{lo^) = ^NM-i{Lj^)) a^, M = 0,...,iV-l 

(8) 

Here the are Gaussian random numbers with zero mean and with correlation 

{a^a^) = S^+^,Q ■ (9) 

This type of 5-anticorrelated noise can be generated rather easily if the symmetry properties of real periodic series 
in the Fourier space (a^) are used 0, Owing to the periodicity of 7(w/i), the index /j, can be made to run from 
—N/2 to N/2, and the maximum frequency is ijJmax = 7r/At p7[ |. The Fourier components are related by a^+pw = 
for any integer p, and = a* because the original random numbers in time are real. The anticorrelated random 
numbers can then be constructed as = -\- ib^ with 6o = 0, Oq = 1, and and 6^ are each Gaussian random 
numbers with zero mean and a variance of 1/2 for /i 0. 

The discrete inverse transform of any sequence ri{oJ^) is then numerically calculated by a standard Fast Fourier 
Transform algorithm. The result is a string of N numbers, ri[ti) which, by construction, have the proposed time 
correlation (Q). However, due to the symmetries of the Fourier transform only N/2 of these values are actually 
independent and the remaining numbers are periodically correlated to them. Thus we have constructed a true 
random process from i = to a maximum time Tcutoff — AtN/2. 

This algorithm is sufficiently general to allow generation of noise of any given correlation function. For instance, 
in we considered Gaussian noise with a Gaussian correlation function, 

lit) = ^^e-* (10) 
TV Stt 

whose Fourier transform can be obtained analytically, 

j{lu) = 2De-^^'^'/\ (11) 
According to our prescription, if we generate a discrete field of random numbers according to 



T-2 \ 1 1/2 



2DNAtcxp \^j^{cos{2nfj,/N) - 1) 



(12) 



for sufficiently large N the resulting correlation function upon Fourier inversion should match Eq. ([l0| ) (it does). 
Similarly, for exponentially correlated Gaussian noise, cf. Eq.(|^), whose Fourier transform is 

- TT^' ^^^^ 

we generate the discrete random numbers according to 

~, ^ / 2DNM \ , , 

"'""'Hi + 

and again reproduce the correlation function (^ upon Fourier inversion. In both cases 2D is the intensity of the noise 
and T its correlation time. We find that it does not matter whether we discretize the oj using the function cos{2TTfi/N) 
as in (^ or sm{iTn/N) as in (|l^. 

A more difficult noise to generate numerically is one with a memory of inverse power law form, 

Here to is an adjustable (small) parameter. For reasons described below and related to our particular algorithm, the 
value of to turns out to depend on /3. When the power law decay exponent (3 > 2 then not only is the intensity 
of the fluctuations finite (and given by 2eto/{(3 — 1), which increases with decreasing /3), but they also have a finite 
correlation time (which also grows as (3 decreases. The effects of this sort of noise are essentially the same as those of 
any Gaussian noise with a finite correlation time, although direct comparison with most of the literature would have 
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to be done carefully because correlation time effects have usually been studied in the context of noise of a fixed 
intensity. When 1 < /? < 2 then the fluctuations have a finite intensity but an infinite correlation time. We have not 
analyzed this case. Even more persistent correlations, and the ones of interest to us here, occur when < [3 < 1. 

In order to implement the spectral method, we would first have to know the Fourier transform of the correlation 
function. The form ( |l5| ) does not have an closed- form Fourier transform and so instead we postulate or guess 
the form 



e/37rw. 



/3 

max 



[it M^^^/N) 



l(l-/3) 



(16) 



where iVmax was defined earlier as n/At and where luq is a low frequency cutoff that is chosen to control the low- 
frequency behavior and avoids a zero-frequency divergence. To check whether Eq. is indeed an appropriate choice 
one has to perform the numerical inverse transform and compare with (p^. For a given exponent /3 this comparison 
involves the three parameters At (the discretization time step), loq, and Iq. The value of loq that leads to the best 
agreement turns out to depend on the time step. In our calculations we have mostly used At = 0.01 but in our 
barrier-crossing calculations in Sec. ^ we take At — 0.02. For this range of At we find that the best choice (giving 
the best concordance between the desired correlation function and the assumed Fourier transform) is loq = 0.0001. 
This is the value used throughout this analysis. 
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FIG. 3: 7(0) vs /3 obtained from the simulations. 
From Eq.(|l|) we obtain an important relation, 

7(0) ^ - / 7(a;)(iw = e + O{ujo/ujmax), 



(17) 



consistent with Eq. (p^. Therefore 7(0) is independent of (3 and is a direct measure of the noise control parameter 
e. This provides a first check for our algorithm. 

A second test of our algorithm is presented in Figs. |^ and where we exhibit the numerically generated correlation 
functions as well as the analytic counterparts with /3 — 0.3 and 0.6. Within the resolution of these figures, the 
simulated results agree well with the expected power laws given by Eq.(^) for the particular parameter choices 
shown in the caption. The value of e has been assumed to be selectable independently of (3 and according to Eq. (|l^). 
Figure ^ shows that this is indeed the case, i.e., that 7(0) as generated from our simulations is essentially independent 
of (3. The small variations are due to the discreteness of the time variable and the stochasticity of the problem. More 
complicated is the choice of the parameter t^: for a fixed e, the choice of t^ that leads to a best match between the 
numerically generated correlation function and the analytic form (|l^) depends on /3, a dependence that is exhibited 
in Fig. 1^. The parameter also varies with the time step At, a dependence we have not exhibited because we hold 
At fixed in our analysis. In any case, implementation of our algorithm requires attention to these dependences. The 
values of t^ indicated in Figs. ^ and |l| result from this best fit. 
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FIG. 4: Parameter to found from the best fitting of the autocorrelation function Eq. (hi 
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0.25 


1.143 


1.151 


0.5 


1.333 


1.344 


0.75 


1.600 


1.596 



TABLE I: Comparison of theoretical and simulation exponents 



III. SUPERDIFFUSION OF A FREE PARTICLE 



The first problem we study is of interest in the analysis of the mean exit time of a free Brownian particle from a 
domain bounded by absorbing walls. We consider a free particle that moves in one dimension under the influence of 
a long range correlated noise until it covers, on average, the distance to one of the absorbing barriers. 

A dynamical equation that describes this type of dynamics is the very simple stochastic differential equation 
(Langevin equation) 



dx , . 



(18) 



where ri{t) is a Gaussian power law-correlated random noise as described in section Note that here and in all 
subsequent examples we "start the clock" at t = 0, that is, in addition to specifying an initial value or distribution of 
X a.t t — we also assume that the noise is "turned on" at that time. Otherwise one would need to be concerned about 
the role of the time t = in the evolution of the noise, which in turn would affect the description of its correlation 
function. 

We are particularly interested in the dependence of the mean arrival (absorption) time on the parameters /3 and e 
and on the length 2L of the interval. The essential scaling features (exponents) of the dependences on these parameters 
can be obtained from a formal integration of Eq. (|l8|) , which immediately leads to the mean squared displacement 
relation 



(19) 



characteristic of superdiffusive behavior . Superdiffusion here arises because the arrival at the boundaries is more 
likely to occur via essentially ballistic motion than in the white noise case, where frequent reversals of dx/dt cause the 
net displacement to be slow. The time at which the mean squared displacement is proportional to can be obtained 
by inverting this expression. The difference between this time and the mean time (T) for first arrival at an absorbing 
boundary is due to recrossings of the boundary and only affects the prefactors, not the exponents. That is, we can 
use Eq. (jl^) to write 

-1/(2-/3) 2^2/(2-/3)^ 



(T) 



(20) 
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FIG. 5: Mean first passage time, (T), as function of tlie interval lengtli, L, for a free particle driven by power law correlated 
noises with different exponents /3. In all cases e — 0.05. The dashed line is the theoretical prediction for white noise. The lines 
joining the simulation symbols are simply guides for the eye. Circles: f3 — 0.25; Squares: P — 0.50; Diamonds: (3 = 0.75. The 
exponents in Table I are fits to these numerical results. 
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FIG. 6: Variation of the mean first passage time with e for a free particle on an interval of length L — 3 and two different 
power laws. Circles: l3 = 0.3; Diamonds: /3 = 0.6. The exponents obtained from fitting the simulation results are 0.61 and 
0.72. The auxiliary lines correspond to the theoretical exponents (1/(2 — /?)), which are 0.588 and 0.714, respectively. 

Note that the (Gaussian) statistics of the noise play no role in this result, which is determined entirely by the 
correlation function of the noise. 

In order to test this theoretical prediction we have performed numerical simulations for different values of the interval 
length 2L, the exponent /3, and the coefficient e. Throughout this paper all Langevin equations are numerically 
integrated using the Heun method, which is an extension of a second order Runge-Kutta algorithm for stochastic 
differential equations The results presented below are averages over 2000 realizations. 

A comparison of the theoretical exponent (2/(2 — /?)) in Eq. ( |20| ) and the simulation results is summarized in Table I. 
The agreement is clearly excellent. In Fig. H the average exit time is presented as function of the interval length on a 
log-log plot. The exit time (T) follows the predicted power law behavior with the slopes of Table I. In the figure we 
also show the result for white noise with its characteristic exponent of 2. 

Fig. ^ shows the dependence of the mean first passage time on the noise intensity parameter e. The simulation 
results are in excellent agreement with the predictions of Eq. (pO|). 
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FIG. 7: Mean first passage time out of an unstable state as a function of e for fixed /? — 0.5. The dashed line is the theoretical 
prediction from Eq. (^6|) and the statistical average is over 5000 realizations. 
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FIG. 8: Mean first passage time out of an unstable state as function of the power law exponent /3 for fixed e — 0.0001. The 
dashed line is the theoretical prediction from Eq. (^6|) and the statistical average is over 5000 realizations. 



IV. DECAY OF AN UNSTABLE STATE 



The decay of unstable states is of course easily triggered by any disturbance, and in many physical instances the 
decay is caused by fluctuations. In its simplest rendition this process can be modeled by a particle initially placed at 
the top {x = 0) of an inverted parabolic potential and subject to fluctuations, 

^ = x-x^+7j{t), x(0)-0. (21) 

The decay at the early stages is dominated by the linear term and the fluctuations, so that in this regime the Langevin 
equation can be further simplified to [^o| 

^=x + rjit). (22) 
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Formal integration of this linear equation immediately yields 

x{t)^h{t)e*, h{t)= I dt'e~''f]{t'). (23) 



The mean first passage time to, say, x — ±1 is identified as the time at which {x'^{t)) = 1. Recrossings are even 
less relevant here than in the previous section; larger values of x require consideration of the nonlinear term in the 
evolution equation. If the time of interest is i 3> 1 (a condition that will be seen to be satisfied by the mean first 
passage time) then the upper limit in the integral in Eq. (23) can be extended to infinity, that is, we set 



x\t)^h'e% dt'e-' r]{t'). (24) 

Jo 

Thus the fluctuations h'^ simply play the role of a random initial condition on x'^ (t) . The random variable h has zero 
mean and variance cr^ given by |2l| , 

POO POO 

= {K") = ( / dt'e-'' rj{t') / dt" e-'" r^{t")) = 



oc 



e 



-2t' 7,/ 



dt' I dt"e-' 7(s) = - t(?e*°r(l - /3, to). (25) 



The stationarity of the noise has been used to calculate the variance; r(l — /3, to) is the incomplete Gamma function. 
The mean first passage time to a;^ = 1 for this process for a generic colored noise was calculated in Ref . : 

{T)=-^-ln{a^) + l (26) 



where 7 = 0.57721... is the Euler constant. This result together with ( p5D gives an explicit expression for the mean 
first passage time. Figures and ^ show excellent agreement between the theory and numerical simulations for the 
dependence of (T) on e and on (3. 



V. THE BARRIER CROSSING PROBLEM 



A problem that is more complex than the essentially linear ones posed so far is the "barrier crossing problem." 
Here we consider the diffusion of an overdamped particle in a double-well potential. The mean first passage time of 
interest is the time it takes the particle to go from one of the potential minima to the other when the transition is 
driven by power-law-correlated Gaussian fluctuations. The particle dynamics is modeled by the following Langevin 
equation: 



dx dV{3L 
H ~ df 



V{t), (27) 



where the double-well potential is. 



1 2 , 1 4 



V{x)^~--x' + -x\ (28) 



Our simulations follow the dynamics of each particle starting in the left well at x{t = 0) = —1.0 until it arrives at the 
right well at x{t = T) = 1.0 for the first time. We then calculate the mean first passage time (T), that is, the average 
of T over many (10, 000) realizations. 

Two questions are of interest: (1) If two noises ri{t), one (5-correlated and the other power-law correlated, lead to 
the same average transition time from one well to the another, are there other properties that allow a clear distinction 
between them? (2) For a power-law correlated noise ?/(t), how does the transition time depend on the correlation 
function parameters? We address both of these questions below. 

The numerical barrier crossing problem is fundamentally different from the problems considered in the previous 
section. In those, no matter how small or persistent a realization of the noise r]{t), the process x{t) will eventually 
reach one of the boundaries of interest. How long a simulation must run in order to make sure that all processes 
have reached the boundaries is essentially a matter of insuring that those with the lowest value of 77 as given in the 
discretization scheme arrive there. In the barrier crossing problem, on the other hand, small-?7 realizations will remain 
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FIG. 9: Typical trajectories of a particle in the double well potential (p7|). Top panel: Gaussian white noise with D = 0.0725, 
(T) = 162.38. Bottom: Gaussian power-law noise with f3 = 0.5, e = 1.15, (T) = 162.87. 
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FIG. 10: Power spectra of the trajectories shown in Fig. ^ . Circles: white noise with D — 0.0725 and (T) = 162.38. Diamonds: 
power- law noise with /3 = 0.5, e = 1.15 and (T) = 162.87. The theoretical predictions from Eq. are also included. 



in one well until the value of 77 changes to a sufficiently high value (easy to estimate ) so as to effectively eliminate 
the barrier. The time it takes to effect this change is of course large (infinite on average) for the fluctuations considered 
herein. In our simulations we have carefully chosen run times that insure that all 10, 000 particles in our ensemble 
have made the passage from one well to the other. As we decrease e (thus weakening the noise) and/or increase /3 
(thus in effect also weakening the noise), the simulation time has to be increased accordingly. Note that these remarks 
address the simulation problem. A complete theory of this process would have to deal with the relative magnitudes 
of a number of infinities (noise intensity, noise correlation, and passage over the barrier of the slowest processes). 

Fig. IP shows two typical trajectories of the barrier crossing process for two different Gaussian noises. One is white 
noise, the other is power-law correlated, and the parameters have been chosen so that the mean first passage times for 
both are essentially the same. The effect of the long range correlations is certainly not evident from these trajectories. 
This observation suggests a more detailed study of the dynamical properties. The first passage time distribution and, 
in particular, the mean first passage time for the Gaussian white noise case, are well known analytically (see pl| and 
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FIG. 11: First panehTypical dependence on e of the mean first passage time from one well to the other. The other parameters 
are /3 = 0.5, uoo = 0.0001, and At = 0.02. The dashed line corresponds to A = 1.933. Second panel: Activation parameter A 
as a function of /3 (the value A = 1.933 at /3 = 0.5 is the result of the slope in the first panel. The dashed curve is the best 
numerical fit with a quadratic form, A{P) = 2.42 + 11.5(1.8/3^ - /?). 




FIG. 12: Mean transition time from one well to the other as a function of (3. The parameters used here are e = 1.15, u:o = 0.0001, 
At = 0.02. 

the many original references therein). Approximate theories for exponentially correlated noise also abound pT] |, but 
to our knowledge there are no results for the infinite correlation time case. To deal with this problem we note that 
the trajectories (in either the white or the correlated noise cases) consist of two distinct components characterized by 
two different time scales. One component describes the rapid fiuctuations in each well, and the other describes the 
much slower switching events between the two wells. We implement this observation by writing the trajectory x{t) 
as a sum of these two contributions: 



x{t) = XD{t) + XB{t), 



(29) 



where Xsit) represents the rapid random motion inside each well, and Xo{t) is a dichotomous random process between 
the values x = 1 and x = — 1, with a characteristic time controlled by the mean transition time (T). The dynamics of 
XB(t) can be approximated by expanding Eq. ( [2^ ) around the minimum of one of the wells, e.g. as x{t) — \ + XB{t) 
in the well. Then x b (t) evolves according to the Langevin equation 



dxB 
~~dt 



-2xB + vit), 



(30) 



which evolves on a characteristic time scale r = 0.5 that is much smaller than (T). 

With this decomposition, the power spectrum of x(t), denoted by S{lu), is just the sum of the spectra associated 
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FIG. 13: Transition time probability distributions for white noise with D — 0.0725. The circles are the results of the simulation 
and the continuous line is the functions Eq. (^5|) with (T) obtained from the simulation data. 




with X£i{t) and xsit). These can be calculated explicitly: 



Si^) - + (31) 

The first contribution is the spectrum of a dichotomous process governed by the time scale (T), and the second 
contribution comes from the exact solution of Eq. ( |30| ) with the noise spectrum 7(0;). In the white noise case 
^{lo) ~ 2D, and for the long ranged noise the spectrum is obtained from Eq. ( p^ ) in the continuum limit, 



e(3iTUj, 



max 



It is interesting to note that the main difference between the white and power law spectra appear in the high-frequency 
(short-time) regime. In the white noise case, S{ijj) ^ whereas in the correlated case we observe the more rapid 
decay 5M - w-^^-/^). 

Figure |l^ shows the spectra for the two trajectories shown in Fig. ^ as well as the prediction (^ij) for each. As 
noted earlier, the parameters have been chosen so as to lead to the same mean first passage time (T), but other than 
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FIG. 15: First panel: Characteristic parameter 9 of the stretched exponential distribution as a function of /3. Secon panel: 
characteristic parameter Ta- The values of the other parameters are e = 1.15 and At = 0.02. 
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FIG. 16: Dispersion of the first passage time distribution as a function of [5 (full circles) compared with the proposed analytic 
form Eq. (pq) (empty circles). 



that there are no adjustable parameters. The agreement between our predictions and the numerical results is quite 
good in both cases. The exponent obtained from the numerical data in the high frequency regime for the correlated 
noise case is 2.49, which agrees with the theoretical prediction of 2.5. It is also clear that a distinction between the 
two types of fluctuations is difficult even at this level of detail for a given mean transition time. 

In the above analysis (T) is the adjustable parameter, that is, we have fixed (3 and e so as to yield a particular 
value of (T) . The dependence of this quantity on the model parameters e and [3 is yet to be explored. 

The first panel in Fig. |l^ is a plot of log(T) as a function of 1 /e. The straight line behavior shows that the transition 
time is described by a law of the Kramers form, 

(T) ^ exp (^^^ . (33) 

Fitting the simulation data in the figure with a straight line we obtain A — 1.933, very far from the value A — 2/27 
obtained for exponentially correlated noise |jri|. We stress that the "activation parameter" A here depends on /3, that 
is, a plot exactly like that of the first panel Fig. [l| but for a different value of (3 leads to a different value of A. The 
values oi A YS [3 collected in this way are shown in the second panel in Fig. |ll|. In contrast, for noises with a finite 
(albeit large) correlation time A is essentially a constant independent of the noise intensity and of the correlation 
time. Comparison with those results is therefore not appropriate. 

The dependence of the mean transition time on [3 qualitatively parallels that of the activation parameter but is 
even more complex because it is non-monotonic. The typical behavior is shown in Fig. |l^. Although we have no 
detailed theory to account for this non-monotonicity, a qualitative explanation is possible in terms of two competing 
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effects, one that causes a decrease in (T) with increasing 13 and the other that causes (T) to increase with increasing 
(3. At very small (3 the mean first passage time is dominated by those realizations of very slow passage for which 
Ty(0) is not large enough to cause a transition and a very long wait is involved before ri{t) changes to a value that is 
sufficiently large. With increasing (3 the wait for change decreases, thus leading to faster passage. On the other hand, 
increasing (3 in effect decreases the intensity of the noise. To see this suppose that we calculate an effective intensity 
by integrating the correlation function only up to a maximum time Unax to), the running time of our simulations: 



to 



1-/3 



eta 



to 



1-/3 



l + O 



(1 - /3)to 

^max 



(34) 



The 1 — /3 in the exponent of the large ratio tmax/to dominates the 1 — /3-dependence in the denominator, so that for a 
fixed tmax this decreases rapidly with increasing /?. This lower effective intensity of the noise leads to slower passage. 
We conjecture that the interplay of the two effects leads to the non-monotonic behavior observed in Fig. |l^, although 
the specific value and specific location of the minimum depend on e. 

In order to verify this minimum and to exhibit more information on the statistics of (T), we have also obtained 
numerical results for the probability distribution P{T) of the exit times (of which (T) is the average). In Figs. |l^ 
and 14 we have plotted the numerical histogram of P{T) for white noise and for a power-law correlated noise, again 
chosen so as to correspond to the same means. The scatter of numerical points around the analytic forms discussed 
below (solid lines) is quite large and could be decreased by doing many more realizations for much longer times. 

In the case of white noise the first passage time distribution is known to be exponential jj, 10 , 



P(T) = — exp 



T 



<T> 



(35) 



Our numerical results agree with this result. For the power law case we conjecture a stretched exponential distribution 
for the transition time, 



p(r) 



TsT{i/e) 



cxp 



T 



(36) 



where the prefactor has been chosen to ensure normalization (r(z) is the Gamma function) and where the parameters 
6 and Ts are obtained from the simulation data. The mean first passage time in terms of these parameters is 



(T) = T, 



r(2/g) 
r{i/e) 



,1 



(37) 



We do not display the e-dependence of the parameters because it is the f3 dependence of the first passage time that 
requires further understanding. Our fitting procedure for the case of e = 1.15 and At = 0.02 gives the /3-dependences 
observed in Fig. |l5|. Although the dependences of the parameters on /3 are monotonic, the dependence of (T) on 9 is 
not monotonic (note that Ts also varies with 0) and gives rise to the observed minimum. Using the data shown in 
Fig. to construct (T) according to Eq. ( ^7| ) one again obtains the minimum as a function of /3 observed in Fig. |lj. 
Another point worth noting is the approach of the exponent 9 to unity with increasing (3 that is seen in Fig. |l5|. The 
exponent is exactly unity for Gaussian white noise and would also be unity for /? > 2. 

From the proposed distribution (|3^) one can evaluate the second moment of the first passage time distribution. We 
obtain 



(T2) T{1/9)T{3/9) 



r(2/6i)2 



(38) 



For white Gaussian noise this ratio is equal to 2. We see in Fig. |l6|that for correlated noise the ratio is greater than 
2 reflecting the greater width of the stretched exponential distribution, but with increasing (3 the approach toward 2 
is evident. The figure shows the moment ratio (T^)/(r)^ obtained from numerical simulations (full circles) and from 
the analytic expression Eq. (^3). 
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VI. SUMMARY AND CONCLUSIONS 

We have presented a numerical study and derived (in some cases) or conjectured (in others) analytic results for 
the first passage time problem in systems driven by long-range correlated noises. We specifically considered Gaussian 
noise, as opposed to noises with Levy or other long-tailed distributions, or noises defined in terms of long waiting times 
between events, that may also lead to long-range correlations. There is an extensive literature on these non-Gaussian 
problems |Q , but very little on the Gaussian counterpart [^2], |lj, ll^ . We have considered highly correlated noise 
with an inverse power law form for the correlation function, Eq. (pa), with < /? < 1. Not only does this noise 
have no finite correlation time, but the correlation function is not integrable. The parameters of the problem are the 
exponent /? and the noise control parameter e. 

First we discussed a numerical algorithm to generate this type of noise. This is in itself not a trivial problem. Then 
we considered several systems driven by Gaussian power-law correlated noises with a special focus on the problem of 
first passage to a prescribed boundary or set of boundaries. Two of these problems, the arrival of a free particle in 
one dimension at either end of a finite interval, and the decay of an unstable state, admit straightforward analytic 
solution. The dependence of the respective mean first passage times on the noise parameters were calculated explicitly 
and compared favorably with numerical simulations. 

Our final application is more complex: Here we considered the first passage of a particle that evolves in a double well 
potential from one well, over the barrier, to the other well. We specifically considered two issues. One is a comparison 
of the trajectories of such a particle subject to Gaussian white noise and to Gaussian power-law correlated noise, 
with parameters chosen in such a way that the mean first passage times for crossing from one well to the other in 
both cases is the same. It is difficult to find a clear difference that would reveal in an experiment which type of noise 
the particle was subjected to. Only the high frequency decay of the respective power spectra shows some differences. 
The other issue we explored is the dependence of the first passage time distribution on the parameters /3 and e. We 
found that the dependence of the mean first passage time on the noise control parameter is of the usual "activated" 
form, but with an effective activation parameter that depends on (3. We also found that the mean first passage time 
exhibits a minimum as a function of /?, which we explain on the basis of competing mechanisms, one of which leads 
to an increase of the mean first passage time with increasing (3 and the other to a decrease. This non-monotonicity is 
a signature characteristic of these fluctuations. We also found that a stretched exponential form for the first passage 
time distribution describes the numerical results, with a stretched exponent 9 and a characteristic time Tg that both 
depend on /3. This stretched exponential form is consistent with the non-monotonic behavior of the mean first passage 
time, and also gives results in agreement with the simulation outcomes for the second moment of the first passage 
time distribution. The stretched exponential first passage time distribution, and the dependence of the parameters of 
this distribution on the noise parameters e and (3, are features that have no analog in Gaussian noises with a finite 
correlation time. 
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